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We present a general quantum metrology framework to study the simultaneous estimation of multiple 
phases in the presence of noise as a discretized model for phase imaging. This approach can lead to nontrivial 
bounds of the precision for multiphase estimation. Our results show that simultaneous estimation (SE) of 
multiple phases is always better than individual estimation (IE) of each phase even in noisy environment. 
The utility of the bounds of multiple phase estimation for photon loss channels is exemplified explicitly. 
When noise is low, those bounds possess the Heisenberg scale showing quantum-enhanced precision with 
the 0(tf) advantage for SE, where d is the number of phases. However, this 0(d) advantage of SE scheme in 
the variance of the estimation may disappear asymptotically when photon loss becomes significant and then 
only a constant advantage over that of IE scheme demonstrates. Potential application of those results is 
presented. 

A general estimation scheme of multiple parameters can be divided into three stages: the preparation of 
some probes, the interaction of the probes with a system which is determined by the parameter vector 0, 
and measurements of the probes after the interaction. Then ff is estimated from the results of the 
measurements. When the dimension of ^is 1, the case becomes single parameter estimation. If the probes are 

uncorrelated, then the central limit theorem states that the estimation error Tr[Cov{ff)] scales as ly/ Vn, with JV 

being the number of resources (photons, atoms) employed. While in quantum world by correlating the probes 
nonclassically, the estimation error may scale as l/ATin an ideal scenario, which is the ultimate limit of precision 
named as the Heisenberg limit' The enhancement in the estimation precision is the main concern of quantum 
metrology, and a lot of work has been done, both theoretically and experimentally^"^''. 

A quantum enhancement in precision is of great importance in metrology such as for imaging and microscopy. 
Recently, the quantum enhanced imaging making use of point estimation theory is presented based on single 
parameter estimation procedure through the Fisher information approach'''. Since phase imaging is inherently a 
multiple parameter estimation problem, the multiple phase estimation is of interest^*. It is found that for unitary 
evolutions, simultaneous estimation (SE) of multiple phases provides an advantage scaling 0(d) in the variance of 
the estimation over individual estimation (IE) of each phase, where d is number of phases to be estimated. This 
conclusion holds for noiseless processes. However, in a realistic scenario, noise cannot be avoided due to 
decoherence. An investigation of whether this advantage still exists for a general evolution is necessary. 

For noisy processes, it is not known in general if and when the quantum enhancement of precision from 1 ^ Vn 

to 1/Ncan be achieved though general expressions for the uncertainty in the estimation are known. The problem 
is that their calculation involves complex optimization procedures. Fortunately, a general framework is proposed 
recently to obtain attainable and useful lower bound of the quantum Fisher information (QFI) in noisy systems^^. 

In particular, this lower bound captures the main features of the transition from the l/Ntol ^ Vn precisions for 

the cases of noisy channels such as photon loss and dephasing. Those results are for the single parameter 
estimation. 

In this work, we present a general framework for the estimation of multiple phases with noise. We apply this 
framework to study a specific example of the photon loss type noise. Photon loss is a very usual noise type in 
optical systems. We make a conjecture that with only photon loss considered the QFI matrix of the phases can be 
saturated for a certain set of initial probes, which means that we are in principal able to find a measurement M to 
make the Fisher information matrix after measurement equal to the QFI matrix. In this way, the QFI bound 
computed is a tight lower bound of the uncertainty of the estimation. We show that in the limit of noiseless, the 
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precision can achieve the Heisenberg hmit 1/N with an advantage of 
0{d) for multiple phase, thus recover the known results^"". With noise 
increasing, SE is always better than IE, but the 0{d) advantage may 
disappear asymptotically, with photon loss taken as an example. At 
the same time, the precision of estimation decreases to the standard 

quantum limit (SQL) 1 ^ So similar as for single phase, our result 

of multiphase can also capture the main features of the transition 
from Heisenberg limit to standard quantum limit. 

We shall consider a multiple phase estimation model described by 
Fig. 1. In the preparation stage, a probe state is created of the form 

D D 



k=l 



We assume that the amount of resources employed in the estimation 
process is restricted by the photon number N, and Ni^ describes the 
kth possible distribution of N photons in different modes, which is 
represented by a vector iNkQ,Ni^i, ■ ■ ■ ,-Nyc where Nk,i stands for the 

number of photons employed in the rth mode ''nd ^ N^: , =N. D 
= {N + dy./N\d\ stands for the total number of possible distributions. 
Normalization is required such that Y^k-i l"^*:!^ = 1- In an estima- 
tion scheme, the probe state is chosen beforehand, and one aim of 
metrology is to find out the optimal probe to estimate the parameters. 
For simplicity we only choose pure states as probes, so we have po ~ 

In the evolution stage, we consider the case that states in different 
modes evolve independently. In the mode /, evolution is determined 

by the parameter Uj, expressed in terms of Kraus operators n[''(f,), 

which satisfies nf ^(^^Af (0,) =11. The evolved state is then 

given by 



where we denote ff 



(Oi, 



nr(x)n;;'(()i)® . 



1 



(2) 



/rf) and 



Results 

The advantage of simultaneous estimation. As is shown^', SE 
provides an 0{d) advantage over IE, without noise considered. 
Here we shall show that even under general evolution, SE is still 
better than IE, but the 0(d) advantage may disappear gradually, 
with photon loss taken as an example. We remark that our results 
of noisy processes can recover the case of noiseless in a continuum 
manner thus possess the SE advantage. 

In Fig. 1, only one reference mode 0 is implemented to estimate the 
d phases 0i to 6^. We now consider the scheme to implement d 
reference modes, with each connected to a corresponding phase. 
The initial state can be written as 



Nk,d) 



(3) 



where each reference mode experiences the same evolution as the 
original mode 0. We remark that any IE strategy is equivalent to use 
an initial state with the form 



llAo) =|l/'l)oi,l®l'/'2)o2,2® ■ ■ ■ 



(4) 



and only separate measurement for each phase is allowed. Now we 
see that IE is actually contained in the complete set of SE strategies, 
which leads to the conclusion that SE is generally better than IE even 
under noise. 

Phase estimation under photon loss. A beam splitter is generally 
used to model photon loss. A possible set of Kraus operators in each 
mode is given by^' 
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Figure 1 | A multiple phase estimation model. An initially prepared probe 
state I {j/o) undergoes a general evolution described by d + 1 sets of Kraus 
operators, depending on d parameters which we are supposed to estimate 
simultaneously. Different modes evolve independently. 



(5) 



where i] is the square of the transmissivity r (ranging from (J = 1, 
lossless case, to i-j = 0, complete loss). It is conjectured in 
Supplementary Material that, as long as all the a*- in Eq.(l) are real, 
for this particular set of Kraus operators, the QFI bound can be 
saturated. Since equivalent sets of Kraus operators lead to the same 
evolved state, the QFI matrix should be the same no matter what 
Kraus operators are chosen. Consider the following set of Kraus 
operators 



(1 -^,-)'' 



(6) 



where 5j are arbitrary real numbers that we are free to choose. 

In the methods part, we have derived a method to give a lower 
bound for the optimal precision of multiple phase estimation 



where the element of the matrix of Cq is 

c„(.,n,)„^4{(^..>.>^-(A<'.)/A«.)J, 

with (■ ■ ■)q standing for s(i/'o| ■ ■ ■ |i/'o)s ^"'^ 

dfl'"'^ 



ll 



= ) 2^1, do, do, ' ' - 



(8) 



(9) 



(10) 



Under the noise of photon loss, following the same calculation as in 
the single phase case, we have^^ 



ajhj + bihi, i=j. 



(11) 



witha; = 1 — (1 + ^,)(1 — !//),&,- = (1 + SjYriiil — fj,). For simplicity 
of calculation, we suppose that t], = rj for all i, or all modes are 
symmetric. 

We first consider the best IE strategy to estimate d phases with 
limited resources of N photons. Generally the minimum uncertainty 
of the estimate of phase can be written as 
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where t is the scahng coefficient under certain conditions with t = 2 
being the Heisenberg scale and t = I being the SQL scale. Q is a 
constant and n, is the number of photons employed in the estimation 
of phase i. Since all modes are symmetric, we assume that under the 
best IE strategy, the uncertainty of each phase follows the same 

scaling. We then need to minimize AO^ = Ct — ■ 

Through basic calculation we know that the minimum is obtained 
when the estimation of each phase uses the same amount of 
resources, which is N/d photons, for any positive t. Then we have 

d ^ 
minVAQf= ' , d. (13) 

Now we turn to the SE strategy. If we choose Sj = 1 and 



1-'?, 



substitute them into Eq.(ll), all the off-diagonal terms of Cq will 
disappear, we then have 



4iy 



-'7(«.>o 



(14) 



from which we can clearly observe the disappearance of the 
Heisenberg scale as is expected. 

To see that the 0{d) advantage may disappear in the asymptotic 
case, we first assume that N^l, d2>l,N/d^l. From Eq.(14), we are 

to seek a state | i/zq) which maximizes . . Since . (m,)q < N, 

1 d^ W 
we have . ^ < — and the equality is attained when («,)q = 



N 

for any i. Then a lower bound for SE is obtained: 
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Figure 2 | A comparison of SE and IE strategies for multiple phase 
estimation with d=2,ffi = 2,ff2 — 2- For (a), r] is fixed at 0.9 and Nis 
various. For (b), JVis fixed at 6 and t] is various. The black solid line gives 
the total variance I AOsEideai^ ^ without any noise using the probe states I ij/). 
The red dashed Kne gives the total variance I A^seI^ under photon loss 
using the probe states I tj/s). The blue dotted line gives a lower bound of the 
total variance lAOj^l^ under photon loss using IE strategy with the optimal 
probe. 



We know in the asymptotic case, the scaling coefficient tin Eq.(13) is 
1, and the total variance is 



d' 



(16) 



Compare Eq.(15) and Eq.(16), we see that the 0{d) advantage no 
longer exists. 

In order to exhibit more clearly the transition from the Heisenberg 
scale with the 0(d) advantage to the SQL scale without the 0{d) 
advantage, we investigate the SE strategy using a specific probe state 
I lAs)- I 'As) is a generalized NOON state as defined in Ref 26, explicitly 
written as 



where a. = 



li/^j) = a(|O,N,0, ■ ■ ■ ,0) + |0,0,N, ■ ■ ■ ,0) + 
|0,0,0, ■■■,N)) + /;|N,0,0,---,0), 

1 



(17) 



and doc^ + [V = 1. The reason we choose this 



state is that in the noiseless case estimation with this state has both 
the Heisenberg scale and the 0(d) advantage^*, and we wUl show how 
they disappear as noise becomes significant. To further simplify the 
calculation, we assume that Sj = 8, which is reasonable since all 
modes are symmetric. Then only one variable S needs to be opti- 



mized to make the lower 
Asymptotically we have 

1 

'4" 



bound Cq as tight as possible. 



rrlc^i 



when 5 — 



N/il 



1-') 



N+l 



1, see Supplementary Material for details. For 



^ — , we have 

17 N 



Tr Cn' = 



1 1 

4{N/df' 

We see that it is the Heisenberg scale, additionally, compared with 
Eq.(13), the 0(d) advantage of SE exists. Whereas for N3> : 



1- 



have Tr I Cq M = 



l-i] 1 



1 



■d. We see that it is the SQL scale and 



4f? N/d 

compared with Eq.(13), the 0(d) advantage of SE disappears. 

Although we have proven that SE provides at most a constant 
increase of precision over IE asymptotically for large noise, it doesn't 
mean that there is no need to use the SE strategy. Rather contrarUy, 
it is shown in Fig. [2] that for d = 2, 1] = 0.9 and small numbers 
N, a significant decrease of uncertainty about 50% can be 
achieved. For IE, an optimized state over all states of the form 

^_^a„\n,N /d — n) is chosen as the probe to estimate an indi- 
vidual phase. We have calculated a lower bound of the QFF^. For 
SE without loss, the state | ij/^) is chosen as the probe. For SE with loss, 
we use the same probe and calculate its QFI matrix numerically. 
Since we have proven that for this initial state, the QFI matrix can 
be locally saturated, we have jA^s^j^ = Tr[lQ(ff)^^]. So in principal, 
an advantage of SE over IE larger than that shown in Fig. [2] can be 
obtained. From the result, we see that if we need to estimate multiple 
phases, we should estimate simultaneously to achieve higher 
precision. 

In Fig. 2, we have also made a comparison of different estimation 
strategies versus various rj. We see that under low ij, which means the 
photon loss is significant, SE using states | xj/s) is worse than IE. This is 
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understandable, because for calculating jA^j^]^, we have used an 
optimal probe, but for calculating jA^<;£p, only is used. Ii/^j) is 
a generalized NOON state and is vulnerable to photon loss. A state 
robust against photon loss may be necessary^'. For higher t], li/zj) is 
enough to beat the IE strategy. 

Discussion 

We have presented a lower bound for the error in multi-parameter 
estimation under noise, within the framework of quantum met- 
rology, and photon loss is exemplified. We have proved the useful- 
ness of this bound by showing that it can capture the main feature of 
the transition from the Heisenberg limit with the 0(d) advantage to 
the SQL limit without the 0{d) advantage as noise becomes signifi- 
cant. We have also shown the advantage of SE over IE in precision. 
The enhancement in precision can also be applied for single phase by 
replicating it to several copies. This novel scheme is better than 
simply duplicating the measurement instrument. Our analysis of 
multiple phase estimation should be of wide interest in many pro- 
blems. Quantum enhanced phase imaging is one potential applica- 
tion. A recent investigation of quantum phase imaging used point 
estimation with single parameter'^, since phase imaging is inherently 
a multiparameter estimation problem, our results provide an 
approach to this problem. Our results should also be of interest in 
gravitational wave detection'^, since it can be recast as optical phase 
estimation'"". They will also motivate an investigation into the role of 
noise in quantum enhancement. Thus, the application of our results 
is worth investigating for various quantum metrology problems. 

Methods 

It is known tliat, the precision of the estimate of ^, described by its covariance matrix 
Cov{ff), is limited by the quantum Cramer-Rao (QCR) inequality'^**'^^ 

Cov{0)>{MlQ{d))-\ (19) 

where the mequality means that Cov(^ — (M7q(^)~^ is positive semidefmite, /q{^ is 
the QFI matrix, M is the repetition of the whole estimation process. Here we have 
assumed that the estimator of ^is unbiased. This is a reasonable assumption since 
Cramer has proved that the maximum likelihood method will give an asymptotic 
unbiased estimate as M ^ ^c^* a brief introduction about the QFI approach for 
quantum metrology is presented in Supplementary Material. Since we are interested 
only in the quantum enhancement, we shall set M to 1 for this letter. The total 
variance of all the phases is then 

d 

\Ae\^ ^ ^ _ Tr[Cov{6)] > Tr[l{0)~^] . (20) 
( = 1 

Inspired by the work^^, we propose a general method to derive an upper bound 
Cq{0Ai) of /q(^, where fli is any Kraus representation of the quantum channel. 
Suppose the real value of the parameter vector is ff, and e is an infinitesimal increment, 
then we have the relation between the Bures fidelity and the QFI matrbc at ff^^: 

{FMO)AO+e)]f = 1 - \Y.<'ieM0hj, (21) 



where the Bures fidelity is defined as: Fb[p,(j] — Tr\J y/p<7y^. Uhlmann's theorem 
states that^** 

{FBlp{e)A0 + ^}]f- max |(a)(^)|T(6» + e))|', (22) 

where |0(^) is an arbitrary purification of in an enlarged space SE, and 
|4:'(^ + e)) runs over all purifications of p(0 + e). Since 

|(<I)(0)|>I'(0-He))|^ =FB(|<t(e)), |«F(e-H€)))^ = 1 - i^__.eie,CQ(e),^, where Cq(^ 

is the QFI matrix at ^in space SE, we have Iq{ 0) < Cq( 0) . The equality may actually be 
achieved. Because for pure states | T(^), its QFI matrix can be explicitly written out. 
This will provide us a method to derive useful analytical bounds of Iq{,0). 

Notice that for the scheme of Fig. 1 , although the probe state may be correlated, the 
evolution is separated for different modes. Thus rather than to purif)' the system 5 on 
the whole, we may purify each mode independently, which greatly reduces the dif- 
ficulty of purification. Add an environment to the respect system S„ and purify^ the 

evolution |n|'' | to a unitary one O^^'^K the evolved state psi^) becomes a pure state 

|T(^)s£, given by 



mO))sE-^^''HdMo)sWE^ (23) 

where tj^^^^O) - U^^^''^ ®'l^^ijf'^'\9i), \Q)^ - ®f=o|0>£,.- The purified unitary 
evolution is connected to the original Kraus representation through the equation^", 

hf (6i) - (/,■ I Uj'^^''^ (0i) |0)£i, (24) 

where 1/,)^ form a basis for the environment 

We show in Supplemental Material that the QFI matrix for the enlarged total 
system SE can then be expressed as 

Go(.,n,),^4{(B^«)^-(A<-.)^(A(..)j, (25) 
with (■ ■ ■)q standing for s{i/'o| ■ ■ ■ |'Ao)s ^^id 

^''-E^'^<' (26) 



^(y) ^ ) l^u d(h dOi ' ^ —J (27) 

So at first place, we have Iq{0) — min^j^ 0^(^,11/), with the minimization running 
over all possible Kraus representations of the quantum channel. In order to reduce the 
difficulty of the optimization process, we only consider independent purification of 
each mode, such that flj(^) -flj^ (x)!!,^ (^)l)(x) • • ■ (g)n,^(fJj). Further we can restrict 
the minimization process to a subclass of all the possible flj, depending on a few 
variational parameters which shall be optimized. The subclass may be constructed 
based on physical insight. In this way nontrivial bound can also be obtained as we will 
present below. 
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